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Abstract 

We extend the dipole formalism of Catani and Seymour to QCD processes involv- 
ing heavy fermions. We give the appropriate subtraction terms together with their 
integrated counterpart. All calculations are done within dimensional regulariza- 
tion. Our formulae can be used with three variants of dimensional regularization 
(conventional dimensional regularization, the 't Hooft-Veltman scheme or the four- 
dimensional scheme). 



1 Introduction 



Studies of hard-scattering processes in QCD with heavy particles are important at present 
or future colhders. Rehable theoretical prediction require at least the evaluation of next- 
to-leading order (NLO) QCD corrections. Next-to-leading order calculations have to 
combine virtual one-loop corrections with the real emission contributions from unresolved 
partons. Treated separately, each of the two parts gives an infrared divergent contri- 
bution. Only the sum of the virtual corrections and the real emission contributions is 
infrared finite. Setting up a numerical general-purpose NLO Monte-Carlo program re- 
quires therefore the analytical cancellation of infrared singularities before any numerical 
integration can be done. The two main methods to handle this task are phase space 
slicing and the subtraction method 0, || . For massless partons both methods are 

available and have been applied to a variety of specific processes, see for example ref. 
and the references therein. For electron-positron annihilation a general formula- 



tion of phase space slicing has been given by Giele and Glover W^. It has been extended 



to initial-state partons by these authors and Kosower [|T^. The extension of phase space 
slicing to massive partons and identified hadrons has been given by Laenen and Keller 



There are two general formulations of the subtraction method. One is the residue ap- 
proach by Frixione, Kunzst and Signer [|T^ and the other the dipole formalism by Catani 
and Seymour [^, both variants can handle massless partons in the final and/or initial 
state. The subtraction method has already been applied to some specific processes with 
massive partons ||16|-|T^. Up to now there is no extension of the dipole formalism to mas- 
sive partons. Dittmaier has considered the dipole formalism for photon radiation from 
fermions [0. In that work infrared divergences are regularized with small masses (as it 
is popular in electroweak physics). However this does not allow a simple application to 
QCD, where divergences are usually regularized by dimensional regularization. 



In this paper we extend the dipole formalism to heavy fermions. The formulae we 
provide are relevant to top, bottom and charm production. With a simple change of the 
colour factors they can be used as well for gluino production in SUSY QCD. Our results 
are not applicable for processes with identified hadrons in the final state, massive initial- 
state partons and processes with different species of massive fermions of unequal masses. 

One-loop amplitudes may be calculated in different variants of dimensional regular- 
ization, such as conventional dimensional regularization (CDR), the 't Hooft - Veltman 
scheme (HV) or the four- dimensional scheme (FD), whereas the Born amplitudes entering 
the real emission part and/or parton densities are given in another variant. Of course, the 
final result has to be scheme-independent. The obvious way to ensure this is to calculate 
everything in the same scheme. Often this is not the most economical solution. Loop am- 
plitudes are most easily computed in the FD scheme, whereas parton densities are given 
in the conventional scheme. Within the dipole formalism it is possible to perform different 
parts of the calculation with different variants of dimensional regularization and to correct 
for the mismatch by some universal terms . We would like to keep this freedom also in 
the extension to massive partons. We therefore keep track of all scheme-dependent terms 
and our final results can be used with any variant of dimensional regularization (CDR, 
HV or FD). 
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Recently, Catani, Dittmaier and Trocsanyi considered the singular behaviour of one- 
loop QCD and supersymmetric QCD amplitudes with massive partons in the dipole frame- 
work [^]. As a byproduct, we confirm their results for the singular terms. 

This paper is organized as follows. In the next section we review the dipole formalism 
and the factorization of QCD amplitudes in the soft and collinear limit. In section 3 we 
outline our calculational technique. In section 4 we derive the D-dimensional dipole phase 
space measure. Section 5 gives the dipole terms together with the integrated counterpart 
if all particles are in the final state. In section 6 we consider the case if there are QCD 
partons in the initial state. Finally, section 7 contains our conclusions. 



2 The dipole formalism and factorization in singular 
limits 

In this section we briefly review the dipole formalism and the factorization properties 
of QCD amplitudes in the soft and collinear limit. We use the notation of Catani and 
Seymour 



The dipole formalism is based on the subtraction method. We explain it for electron- 
positron annihilation, where all QCD partons are in the final state. The NLO cross 
section is rewritten as 



^NLO ^ I j^R , / j„V 



n+1 n 

= J (da^^ - da^) + J Ida^ + J da' 

n+1 n \ 1 

In the second line an approximation term da^ has been added and subtracted. The 
approximation dcr^ has to fulfill the following requirements: 

• da^ must be a proper approximation of da^ such as to have the same pointwise 
singular behaviour (in D dimensions) as da^ itself. Thus, da^ acts as a local 
counterterm for da^ and one can safely perform the limit e 0. This defines a 
cross-section contribution 



NLO 
'^{n+l} 



{da"" 



- da^ 

e=0 



e=0 



n+1 



• Analytic integrability (in D dimensions) over the one-parton subspace leading to 
soft and collinear divergences. This gives the contribution 

<f = j(da^ + JdaA . (3) 

n \ 1 / 

The final structure of an NLO calculation is then 

„NLO „NLO I „NLO ( a\ 

(J = cr{„+i} +^{n} ■ (4) 
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Since both contributions on the r.h.s of eq.(^ are now finite, they can be evaluated with 
numerical methods. The (n+ 1) matrix element is approximated by a sum of dipole terms 



pairs i,j kj^ij 



E E \ fa')' ^' ■■■l ^'^rr.^'^ ViJ,k\'^^ fa'), ■■; k, ■■■), 

pairs i,j k^i,j "i «i 



(5) 



where the emitter parton is denoted by ij and the spectator by k. Here Tj denotes the 
colour charge operator |T^ for parton i and Vij^k is a matrix in the helicity space of the 
emitter with the correct soft and collinear behaviour. \1, (ij), k, ...) is a vector in 
colour- and helicity space. By subtracting from the real emission part the fake contribution 
we obtain 



da^-da^ = d(f)n+l[\M{pi,...,Pn+l)\'^e'^_l^{pu...,Pn+l] 



- J2 H^ihk{pi,-,Pn+l)On\Pl,-,Pij,-,Pk,-,Pn+l] 
pairs i,j k^i,j 



(6) 



Both da^ and da"^ are integrated over the same (n + 1) parton phase space, but it should 
be noted that da^ is proportional to whereas da"^ is proportional to Here 
denotes the jet-defining function for n-partons. 

The subtraction term can be integrated over the one-parton phase space to yield the 
term 

I(^da^ = J da'^= J2 J d<PdipoieVij,k- (7) 

^ pairs i,j kj^ij 

The universal factor I still contains colour correlations, but does not depend on the unre- 
solved parton j. The term I (g) da^ lives on the phase space of the ra-parton configuration 
and has the appropriate singularity structure to cancel the infrared divergences coming 
from the one-loop amplitude. Therefore 

da^ + 1® da^ (8) 

is infrared finite and can easily be integrated by Monte Carlo methods. 



In order to extend the dipole formalism to massive fermions, we have to provide three 
ingredients. First, we need the correct subtraction terms for the real emission contri- 
bution. Second, we have to integrate these subtraction terms over the dipole phase 
space. The integrated terms are combined with the virtual corrections. Third, we have 
to specify a mapping of the momenta, which relates the (n + l)-parton configuration 
Pi,...,Pi,...,Pj,...pk,...Pn+i to the n-parton configuration pi, ...,pij,...,pk,...,Pn+i- 

In order to find appropriate dipole terms, one considers the soft and collinear limits 
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of the matrix element. In the soft limit where the momentum of parton j becomes soft, 
the m + 1-parton matrix element behaves as 



(1, ...,m+ 1|1, ...,m + 1) = -An/jL^^asil, m| J'^^J^]!, ...,m) 



(9) 



with 



T'^'^T II 



m) 



2E 



Am: 



4E 



( '2-PiPk 



T.-T»+EpJ^T..T.||l,...,™) 



We used 



^ V ) (2PiPi + '^PjPk) C^PiPj ) ^ 



Ti •Tfe|l,...,m). 



{2piPj){2pjPk) {2piPj){2piPj + 2pjPfe) (2pjPfe) (2pjPj + 2pjPfc) 
and colour conservation 

^T,|l,...,m) = 0. 



(10) 
(11) 
(12) 



The colour charge operators Tj for a quark, gluon and antiquark in the final state are 

quark : {...qi...\T^^-\...qj...), 

gluon: {...g^...\ir%.g\..), 

antiquark : {...qi...\ (-Tfj) \...qj...). (13) 

If the particles are massless, there is also a singularity in the coUinear limit. For final-state 
partons the momenta are parameterized as 



Pi 



Pj 



zp + k±- 



n 



z 2pn' 



1^1 - z)p -k±- 



1 — z 2pn 

Here n is a massless four-vector and the transverse component k± satisfies 2pk± = 2nk± 
0. The coUinear limits occurs for k^, — > 0. In this limit the matrix element behaves as 



(14) 



(1, ...,m + 1|1, ...,m + 1) = A-K/i^^asil, ...,m\--^P(^ij)^i{z,k±,s)\l, ...,m 



The splitting function are given by 

{s\Pqq{z,k±,e)\s') = 5ss'Cf 

{s\Pqg{z,k^,e)\s') = Sss'Cp 
{l^\PUz,k^,e)\iy) = Tr 
{li\Pgg{z,k±,e)\iy) = 2Ca -g 



2z 



1-z 
2(1 -z 



'^PiPj 



+ {l-pe){l-z) 



(15) 



+ (1 - pe)z 



g^'' + Az{l - z 



1-Z 
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+ 



kl 
1-z 

z 



2{l-pe)z{l-z[ 



.(16) 



We introduced the parameter p, which specifies the variant of dimensional regularization: 
p = 1 for the CDR/HV schemes and p = for the FD scheme. Later on we will chose the 
dipolc terms to have the same soft and coUinear behaviour as the appropriate limit of the 
(n + l)-parton matrix element. 

If the emitting particle is in the initial state the coUinear limit is defined as 

Pa = P, 

, X , I Til 

Pi = (l-x)p + A;x-- — , 

\ — X 2pn 

kl Tl / -I r7\ 

Pai = Xp-k_L + . (17) 

1 — X 2pn 

The colour charge operators for a quark, gluon and antiquark in the initial state are 

quark : \...qj...), 
gluon : {...g^...\tr\..g\..), 
antiquark : {...qi...\T^j\...qj...). (18) 

We denoted in the amplitude an incoming quark as an outgoing antiquark and vice versa. 

3 Calculational technique for the integration 

The integration of the dipole terms over the dipole phase space is highly non-trivial. We 
first find a suitable parameterization of the phase space such that all integrals are of one 
of the following types: 

1 

j (ixx^-^l - xf-"-"- -- 



ldxx''-\l-xr-''-\l-xox)-'' = £(^ff|£^2F,(a,6;c,a;o), 

^ ^ 

[dxx''-'{l-xy-''-'{l-x,x)-''{l-X2x)-'' - 3^)I^^Jl^Fi(a,6i,62;c;xi,X2). 

J r(c) 

(19) 





1 



r(a)ri 


[c- 


a) 


r( 


c) 




r(a)ri 


[c- 


a) 


r( 


c) 




r(a)ri 




a) 







These integrals yield the Eulcr-Bcta function, the hypcrgcomctric function and the first 
Appell function, respectively. The last two are then rewritten as a Taylor series in xq or 
xi and X2, respectively. 

2Fi{a,b;c,Xo) = E ZJ^ 

n=o y^jn "'■ 

F^{a,h,b2;C;Xi,X2) = £ £ {a)rn,+rn,{blUAh)rn, xT^X^^^ ^20) 



mi=0 m2=0 



{c}mi+m2 mil ma! 



Here (a)„ = r(a + n)/r(a) denotes the Pochhammer symbol. Taylor expansion allows 
us to perform the next integral, which is of the same type as above. After all integrals 
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have been performed, we end up with multiple sums involving F-functions depending on 
e = {A — D)/2. These F-functions are then expanded according to 

T{n + e) _ 
T{n)V{e) ~ 



e{l + eZ^{n - 1) + e'Zuin - 1) + £^^iii(n - 1) + ... + e'"^Zu...i{n - 1) 



where Zii i{n) are Euler-Zagier sums defined by 



(21) 



Zmi,...,mk{'^) 



V — — — 

ii>i2>...>Jfc>0 '1 '2 'fc 



(22) 



Rearranging summation indices we recognize that all sums fall in the class of Goncharov's 
multiple polylogarithms p3|, 



Li 



mj.,...,m2,ni 



1 (a^fc, X2-, Xi) 



^1 ^2 

■mi ■m2 ■ ■ ■ -rrik ■ 



n>i2>...>ifc>0 ^1 ^2 

In the case with only one massive parton it is actually sufficient to restrict oneself to 



(23) 



harmonic polylogarithms |2^, defined by 



E 



jl>j2>...>jfe>0 '1 ^2 



1 1 

■nil „-m2 ■ ■ ■ ^-wifc ' 

61 i 



(24) 



The additional dipole corresponding to gluon emission from a massive quark-antiquark 
system can be expressed in terms of two-dimensional harmonic polylogarithms E7\ . 



4 The dipole phase space measure in D dimensions 

In this section we derive the appropriate phase space measure for the dipoles with massive 
particles. Since singularities are regulated by dimensional regularization, this has to be 
done in D dimensions. The phase space measure for n particles in D dimensions is given 

by m 



with 



(27r)^5^ P _ y p, TT \ (25) 



= \l{p,f + ml (26) 



The phase space measure factorizes according to 

d(j)n{P Pl,...,Pn) = -^dcpn-j+liP ^ Q,Pj+l, ■■■,Pn)dQ'^d(f)j{Q ^ pi, ...,Pj). (27) 



4.1 Phase space measure with no initial-state particles 

We first evaluate the two particle phase space 

UP^P^.,Pk) = I ^Snlt'r. .olt-foz. WS^iP - P., - Pk) (28) 



in the rest frame of P, e.g. P = {VP^, 0). We obtain 

where ^d-i parameterizes the solid angle of the {D — 1) spatial components of pk in 
[D — 1) (spatial) dimensions. With this convention we have 

dno = , , (30) 

We then evaluate the three particle phase space 
h{P ^ VuPvPk) 

/di Pi ^ Pi ^ Pk / \ j~) r- j~) / \ / \ 

(2rr)-~^2K(2n)--^2EA2n)o^^2E^ ^''^ ^ - - ' (31) 



in the rest frame of P, e.g. P = (vP2,0). We shall orient our frame such that the 
solid angle of the spatial components of pk coincides with the solid angle of the spatial 
components of pk in (^). It will be convenient to parameterize the spatial components 
of Pi with spherical coordinates, using pk as polar axis, e.g 

d^-^Pi = \pi\^-^d\p,\deism^-^eidQ^S-2- (32) 

We therefore have 

2piPk = 2 \pi\\pk\ cos 9i. (33) 

Finally we obtain the three-particle phase space as the product of a two-particle phase 
space and a dipole phase space: 

03(P Pi,Pj,Pk) = / d(t)2{P pij,Pk)d(j)dipoie{'2piPj, 2pjPk, 2piPk, ^0-2)1 (34) 



where 



d(t)dipoie{'^PiPj, '^PjPk, 2piPk) = {27ry-^-j^ y- N P^ -rri^i - ml^ - ^m\m\ 

{-P''\{^p^,\\p,\\\Pkf)) ' e[-\[%\\\p,\\\Pkf))- (35) 



Note that we already performed the angular integration over dQ^^_2- The triangle function 
A is defined by 

X{x, y, z) — x'^ + y'^ + z'^ — 2xy — 2yz — 2zx. (36) 

We further have 

^ { d2 2,2 2 n \^ 2 

— [P -rrii+mk- rrij - 2piPj) - m^, 



I -' |2 
\Pk\ 



" " ^(^'-^?-^fe + ^?-2p^P/c)'-m5 (37) 



so that 

-p'x (ipiiMp.iMpir) 

= '^PiPk'^PiPj'^PjPk - ml {2pjPkf - ml {2piPjf - m] {2piPkf + ^mlm]ml. (38) 

To find a suitable parametrization of the dipole phase space, we treat the various cases 
separately. We consider first the case where two particles are massless rrii — rrij — 0. 
With D — 4 — 2e the dipole phase space becomes 



d(l)dipoie = T{l-e) ^ ^ ^^"-^ 
1 

• J duu^~^%l - uy-^' (1 - Uou)-^+' J dvv-'{l - v)-', (39) 



1 1 



where 



u 



p2 ' 

2piPj + 2pjPk _ 2piPj + 2pjPk 



2piPk + 2piPj + 2pjPk P'^ -ml 
^ ^ 2piPi {2piPk + ml) ^ 2pjPj {2piPk + m|) 

2pjpfc (2piPj- + 2pjPk) 2piPk (P2 - - 2piPk) ' 

In the case where one mass vanishes {mj = 0) and the other two are equal {mi — m^ — m) 
we obtain 

(47^)'"' /„2\l-^. .2-2e^2e-l 



■'dipole 



r{l-e) 

1 1 
• J drr^-^'{l - r)-^+^/2 _ ^^^)-V2 j {l - s") , (41) 

-1 

with 



2 



Vo( 


;i 


-r) 






ror 




;i 


-r) 






ror 



2p,Pfe = Vror|l + g^/ ':;^\y |. (42) 



4.2 Phase space involving initial state particles 

If initial state particles are involved, we obtain the following convolution: 

1 

d(j){pa + Pb ^ K + Pk + Pi) = j dxd(j){pa+ Pb ^ K + Pk)d(f)dipole (43) 



with 

Pa = XPa- (44) 

In the case = mj = we obtain 

d(l)dipole 



1 



T{l-e 
with P = Pk + Pi - Pa and 



[-P^y "xl-'x~^+%l-xy~^%l-Xox)-^+' J dww-'{l~w)-' (45) 





-P2 



Xq 



X 



171% — 

2paPi + 2paPk - 2piPk 
2paPi + 2paPk 



2paPi{2piPk + ml) 
w = -. (4o) 

2piPk [2paPi + 2paPk) 



Eq. (^) and eq. ( ^Bj ) are derived as follows: From the factorization of phase space we 
have 

d<p{pa + Pb^ K + pk + Pi) = ^d(f){pa + pb-^ K + Q)dmQd(f){Q ^ Pi + Pk). (47) 

Ztx 

We then derive 

d(f){pa + Pb ^ K + Q) = xd(j){pa + Pb ^ K + Pk). (48) 

This equation can be obtained by writing out the explicit expressions for the phase space 
measures in the rest frame of pa + Pb and pa + Pb, respectively. These two frames are 
related by a boost. Further, since pa = xpa, the boost does not affect the transverse 
components. Next, one writes out the parameterization for d(f){Q — > Pi + Pk) in terms of 
sohd angles of particle pk, singles out one angle 9i = l(pk,Pa) and replaces the integrations 
over drnqd cos 6i by integrations over dxdw. 

5 The dipole terms for final-state singularities with 
no initial-state partons 

In this section we give the dipole factors corresponding to the case where all relevant 
partons are in the final state. Initial-state partons are considered in the next section. 
We distinguish the cases of (i) a massless emitter and massive spectator, (ii) a massive 
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emitter and a massless spectator and (iii) a massive emitter and a massive spectator of 
equal mass. We follow closely the notation of Catani and Seymour . 



The generic form of the dipole terms is given by 

1 



-^Pi ■ Pj ^ ij 

The explicit forms of the functions Vij^k are given below. 

5.1 Massless final emitter, massive final spectator 

We consider first the case of massless emitter (particle i) and a massive spectator (particle 
k), both in the final state. The variables y and z of Catani and Seymour are given by 



y 



2piPj 



'^PiPk 



2piPj + 2pjPk + 2piPk 
and related to the variables u and v as follows: 

u{l — u) 

y = uqv 
z{l-y) : 



1 — UqU 

1-u. 



Uq=1 



2piPk + 2pjPk 



ml 



{Pi + Pj+PkY 



(50) 



In the collinear limit eq.(|l^) we have 

y ^0, z ^ z 

and in the soft limit pj — ^ 

y^O, z^l. 



(51) 



(52) 



(53) 



The subtraction term for the splitting g ^ gg has in addition to match the soft limit 
Pi ^ 0, corresponding to y ^ 0, z 0. 



As subtraction terms we use 



2^(1-1/) 
l-z{l-y) 

4 



+ (l-p£)(l-^(l-2/)) 



2piPj 

< z{\-y) 



{\-z){\-y) 



\-z{\-y) l-{\-z){\-y)^ 



where the spin correlation tensor is given by 



(54) 



{z{\ - y)p>: - (1 - z{l - yM){z{l - y)p^ - (1 - zil - yM). (55) 
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The momenta are mapped as follows: 



Pi = api + bpj + cpk, 

Pk = (1 - a)Pi + (1 - b)pj + (1 - c)pk, 



(56) 



with 



yuo 



{2u-l- y{l - u)yuo + 4:u{l - u)y 



2u{l-u) - 



N 



1 

N 



1-u- 



-2yu - uo ((1 - 2m - u'^)y + 1 - 3u + 2^^)] | , 

2y{l -u) - Uo ((1 - u)y'^ + (1 - -u + u'^)y + u- 2m^) 



N = + {1 - uf + {1 - u)y. 
This mapping satisfies momentum conservation 

Pi+Pk = Pi+Pj+Pk, 

and the on-shell conditions 



In the limit y — > we have 



In addition we have 



~2 n ~2 2 

Pi =0, Pk = rrik. 



limpi =pi+ pj, linipfc = p^. 



s,.Pt = 0, 



(57) 



(58) 



(59) 



(60) 



(61) 



e.g. Pi is orthogonal to the spin correlation tensor. The integral over the spin correlation 
can be written as 



j d(j)d 



{2pip,f 

Using S^yp\ = and p^ = this reduces to 

C2ipm - C2, (-9^" + 2 



S^"" = C2iPm + C22ptPl + C2s{PtPl+PtP'^) + C,,g'^''. (62) 



Pi Pk+PkPf 
'^PiPk 



(63) 



Due to gauge invariance only the term —C24:{—g'^'') will give a non-vanishing contribution. 
C24 is obtained by contracting with g^^i,: 



24 



'dipole / n \ 2 diJ.i' 



i^PiPjY 



(64) 
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Integration of the subtraction terms yields: 



2nT{l-e) \{P' - ml) 



2\2 ) 19 



Vgg{U0,e), 



Vgg{U0,e) 



1 1 



= Cf5ss' JduJ dvu-^'{l - u)-^'{l - uouYv-^-'(l - v) 

T{-2e) 







2 + {l-pe)u 



u 



C,5..,E(=£lI(l^r(2 - 2.) . , , 

T{l-2e) ^ ' [ r(2-4£)- 

+ (1 - pg) p|3l^^| 2i"i(-£, 2 - 2£; 3 - As; u^) 

1 3 17 1 5 
£2 + 2£^T^2^~ 6" ' 2uq ' 2 
+2Li2M]+0(£), 

as 1 / lA//2p2 y 



Fi(-£,-2£;2-4£;iio) 



— CpSss' 



2 1 1 (1 - 'Wo)(l - 3mo) , X 







qiqj,k 



27rr(l-£) \{P^-mlf 



^qq{UO,e), 



Vqq{uo,e) 

1 1 







2piPj 



r(-£)r(i-£) 



r(l - 2e) 

2 T{2-2eY 



r{l-2ey 



1 - pe r(4 - 4s) 



r(2 - 4s) 
2Fi(-£,2-2£,4-4£,Mo) 



2Fi(-£,l -2£,2-4£,Mo) 

+ gauge terms 



2 13 1 2(1 -Mo) 
-— - — + 77P + 



3s 6 3 3 Mq 
^ 1 (l-_uo) / 2 _ „^ + 2) ln(l - Mo) 

O Mo 



a„ 1 



+ gauge terms +0{s), 



Vgg{uo,s) 



1 1 



= 2Ca JduJ dvu-'^%l - u)-^%l - UquYv-^-%1 - v) 







ln^ + - ^ ) + (1 - p^)^^/^- 

M 1 — M (1 — m) (1 — MoM(1 — M))y '^PiPj 



r(-e)r(i-s) 
r(i - 2s) 



r(-2s)r(2 - 2s) 

r(2 - 4s) 



(2Fi(-s,-2s,2-4s,Mo) 



T(2 - 2eV 

+2Fi(-s, 2 - 2s, 2 - 4s, Mo)) + ^4 _ ^^-^^ 2 - 2s, 4 - 4s, mo) 

^or ( n^'^\( . ^ ^(l-s)^(-2s) - - r(n + l-s) r(m + n + 2-2s) 
^'""^^"^ r(-s) r(n + 2 - 2s) r(m + n + 2 - 4s) 

13 



r(m-£)^ „ 



ml 



+ gauge terms 



1 11 1 5 2 67 1 (1 - Mo) .... X 

e-^ 6 e 6 12 3 Mq 



(2 -1^0) + 11 In(l-'Uo) 



6 



Mo 



+ gauge terms + 0{e). 



(65) 



5.2 Massive final emitter, massless final spectator 

We now consider the case of a massive emitter (particle k) and a massless spectator 
(particle i) . It is sufficient to consider the case where a heavy quark (or antiquark) emits 
a gluon. We now have 



y 



and 



y 

z{l-y) 



2piPj + 2pjPk + 2piPk ' 



/ (I-m)mo ^ 
It 1 V 

\ 1 - UqU , 



2piPk + 2piPj 



(66) 



Mq = 1 - 



ml 



{Pi+Pj+Pky 



l-u. 

The momenta are mapped as follows: 

emitter : Pk = Pk + Pj 



(67) 



i-y 



Pi 



spectator : p^ 



1 



l-y 



Pi- 



rn 



Since there is no coUinear singularity, we just have to match the part of the soft singularity 
which corresponds to this dipole factor: 



{^\'^<lk9j,i\^ 



2z{l-y) ^fl-uo\l 



With 



l-z(l-y) 

Uo{l - u) 



uq J y 



1 — UqU 

we obtain for the integral over the dipole phase space 



(69) 



(70) 



d(j). 



1 



'dipole 



"^PjPk 



a. 



qk9j, 



27rr(l-£) \(P^-mky 



VQp(Mo,e), (71) 



1 1 





uq{1 -u) - [l- Uq){1 - vqvY'^ 



14 



. 2r(l -e)tt nm-2enn + m + l-e) 
^to^o nn + m + 2-As) 

n + 2-2s . ^..^ , 

-uq - (n + 1)(1 - Mo) 



— CpSss' 

+Oie) 

Here we used 



n + m + 2 — ie 
1 



n+m 
"O 



ml 



(1 + ln(l - uo)) + 4 + ln(l - uo) - 4Li2(Ko) - 7; In^l - Mo) 



00 ii-l h 1 



Lii,i(l, lio) = Hu{uo) = - ln^(l - Uq), 



00 ii — 1 



n=2 12=1 



«2 



Lii,o(l,Mo) 



Mo 



ln(l - Mo). 



(72) 



(73) 



5.3 Massive final emitter and massive final spectator of equal 
mass 

We now consider the case of a massive emitter (particle i) and a massive spectator (particle 
k) of equal mass {rrii = ruk = m). It is sufficient to consider the case where a heavy quark 
(or antiquark) emits a gluon. Since there is no coUinear singularity, we just have to match 
the part of the soft singularity which corresponds to this dipole factor: 



•'\ 



{2piPj + 2pjPk + 2piPfc)2 - {2m'^y 



{2piPkY - (2m 



2^2 



2piPfc 



- 2- 



m 



1 



r)(l — Tor) 



2(l-ror)-(l-ro)- 



2piPj + 2pjPk 2piPj 
1-ro 



sos. 



(74) 



where 



s — 



Sq = 



2piPj + 2pjPk 



2piPj + 2pjPk + 2piPk - 2m^ ' 

1 2piPj - 2pjPk 



So 2piPj + 2pjPk ' 



/ ro(l - r) 
1 — Tor 



, ro = 1 



4m^ 



{Pi + Pj + Pkf 



(75) 



The singularity occurs for r — > 0. In this limit, the expression in the square root tends to 
1. The inclusion of the square root term facilitates the analytic integration of the dipole 
term. The momenta are mapped as follows: 



Pi 
Pk 



P- 
P + 



ro 



yl 



:i-ro) 



^^0 



To) 



(76) 
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with P — Pi + pj + pk and 



yo = p-2 ■ (77) 

This mapping satisfies the on-shcll conditions pj = p\ = fn^ ■ In the soft hmit Uq tends 
to 1 and the correct asymptotic behaviour of the mapping is easily verified. Integration 
yields: 



gi9j,k 



- I '^^<^voiei^V,^9^,k - ^ f^T^ ) ^Qq(^o, el (78) 



VQQ(ro,£) 



{l-r)-'{l-ror)-^ / rfs(l - s^) 



[(2(1 - ror) - (1 - Tq)) (1 - Sos)"' - (1 - ro)(l - Sqs)-' 

1 OC CxD OC 

= c^5..,-^22 Y: u-^y (i + i-iy) ri^^" 

^ 1=0 j=0 ifc=0 

r(l + i + j)V{2 - e) V{i + e) V{k + j/2)V{k - 2e)V{l + j72 - e) 
r(2 + i + J - £) r(£)i! V{j/2)V{k + j72 + l - 3£)A;! 

2-(l-ro)(2 + ,y-±^ 

- C.^..|-fl-^^ln^)-21nro 
{e\ 2 ^ 1 - 



^ + ^° ^Li2 (V^) - Li2 (- V^) + 2 Li2 f- 



2^A^ 



+Li2 



ro - 1 



2v^ 



- Lio 



ro - 1 



2 
1 



2 Li, 



-2 Inro In ( i±v^ ] + In 2 In 



Lig 



- Li2 



1 - y/n)^ 



1 + 2 



+ ln(l - V^) In + \ ln'(l + V^) - ^ ln'(l - V^o) 

+0(£). 



(79) 



Here we used 
I] I] = Lii,i(a;2,a;i) 

= ln(l - xi) ln(l - 0:2) + Li2 



X2 



I — X2 



- Lio 



2^2(1 - xiY 



I — X2 



-- In^ (1 - X1X2) + In (^2(1 - Xi)) ln(l - X1X2) - IA2 



1 - X2 \ 
1 - X1X2 



j + Li2 (1 - X2) , 
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oc 12 
jl=2i2=l 2 

oo n-1 „n 

11=212=1 '1 



•^2 



Lii,o(a;2,a;i) 
Lio,i(a;2,a;i) 



Xi 



1 — Xi 
1 



ln(l - XiX2)-, 



I- X2 



ln(l — X1X2) — 



X2 



1- X2 



ln(l — Xi). 



(80) 



6 The dipole terms with initial-state partons 



In this section we consider initial-state partons. We distinguish the cases of (i) a massless 
emitter in the initial state and massive spectator in the final state, (ii) a massive emitter 
in the final state and a massless spectator in the initial state. The generic form of the 
dipole terms for case (i) is given by 



2pa • Pi X 
In the case (ii) we have 

1 1 



1 1 , ~ , 

-(l,...,(a«),...,/c, ...| 



T T 

T2. 



Vni,...,iat),...,k,...). (81) 



{l,...,{ki),...,d, 



-■- a ■ ■'-ki 

rF2 
^ki 



V:,\l,...,{ki),...,~a,...). 



(82) 



2pi • Pk X 

The explicit forms of the functions V)?* and V)," are given below. 

6.1 Massless initial emitter, massive final spectator 

We now consider the case of a massless emitter (particle a) in the initial state and a 
massive spectator (particle k) in the final state. The variables u and x of Catani and 
Seymour are given by 



u — 



'^PiPa 



2piPa + 2pkPa ' 

The variables u and w are related by 

1 — X 



X — 



2piPa + 2pkPa - 2piPk 
2piPa + 2pkPa 



u 



1 — XqX 

As subtraction terms we use 



W, Xq — 1 — 



mt 



1 — x + u 



ml - {Pk+Pi-Paf 
) 

(1 + x) - pe{l - x) 



(83) 



(84) 



87r/x^^asTR(5ss' [1 — — 2x{l — x)] 



X 



1 — ,T + n 



2piq 

1 + a;(l - x)j 



with 



X 



1 M 1 

-Pi - ( 

r 1 — r 



2piq 



1 . 1 

-Pi - ' 

r 1 — r 



(85) 



(86) 
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Here q is an arbitrary null vector not equal to pi and 



The momenta are mapped as follows: 

Pa = XPa, 

Pk = Pk + Pi - {'^ - x)Pa- 

This mapping satisfies momentum conservation 

Pk-Pa = Pk+Pi- Pa, 

and the on-shell conditions 

Pa = 0> Pk^ ^k- 

The desired asymptotic behaviour 

linipa = xpa, linipjfc = pk 



(87) 



is fulfilled. In addition we have 



« = 0, 



(89) 



(90) 



(91) 



(92) 



e.g. Pa is orthogonal to the spin correlation tensor. We start with the integration of V^^"^'. 
We obtain 



dipole 



1 y,ag. ^ ^-^ 1 



47r/i 



'iPiPa 



27rr(l-£) Vm|-P2^ 



^ J d(t), 

V««(x,xo,£) 

1 

= x%l - x)-^^(l - x^xf j dww-'-^{l - w)- 



2 \ <^ 



V^^(x,a;o,£), 



Cp^ss' 



1 — X -\-u 



[1 + x) - pe{l - x) 



C}<^6ss'X~{l — x) ^(1 — XqX) 



.V{-e)T(l-e) 



1 — X 



r(i - 2s) 

1, -e; 1 - 2e; ~^ ) - {1 + x) - pe{l - x)) . (93) 

1 — XnX J J 



V^"^ is a distribution in x. In order to obtain expressions which are integrable at a; = 1 we 
rewrite 



1 

;i - x)-^^-' = (1 - x)-^^-\ + 5(1 -x)j dy{l - y) 



-2e-l 



(94) 



where 



{l-x) 



-2£-l 



1 — X 



2e 



ln(l — x) 
1 — X 



(95) 
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In order to expand the hypergeometric function we use the analytical continuation formula 



2Fi(a, b; c; -x) 



(x) f(W^^^H"'^ + ""''^ + ""^'"x 

T[a)T[c — b) V x 

largxl < TT. 



(96) 



We obtain 

V''''(x,Xo,£) 



6 \ 1 — X 



(1 + x) +p(l-x) 



+5(1 ~x)(\ + - ln(2 - xo) + ^ + 2 ln(l - xq) ln(2 - Xq) + 2 Li2(xo - 1) 



e 



6 



1 



- ln^(2 - xo) - 2 (Inx + ln(2 - Xqx)) , 
2 / \ 1 — X 



ln(l — x) 



+ / 1 - a; 



(97) 



+ (1 + x) (Inx - 21n(l - x) + ln(l - xqx))] + 0{e). 
For the integration of the other dipole terms we proceed in a similar way and we obtain 

1 / 4V ^' 



d<Pdipolp 7^ V/"''* 



V^«(x,Xo,£) 

~2ei 



V^^(x,xo,£), 



x'il - x^il - xqxY J dww^'~\l - w)-'Tr6, 



Tn6ss'X%l - x)-^%l - xox) 



27rr(l -£) \ml - 

1 — pe — 2x(l — x) 



1 — pe 

r(-£)r(l -e)l- pe- 2x(l - x) 



r(l-2e) 



1 — pe 



TrSss' 



(x^ + (1 - x)^) + 2px(l - x) - (x^ + (1 - x)^) (Inx - 21n(l - x) 



+ ln(l - Xqx))] 
+0{e), 



{9t 



1 ^s{9)^C^q^ _Ots 1 
(^(Pdipolp- —Vi. 



27rr(l-£) \ml-P\ 



V''^{x,Xo,e), 



V''^(x,xo,£) 

1 

= x''(l - x)~^^(l - xqx)^ J dww~^'^{l - w) 



C, 



X IViq 

r(-e)r(i-5) 



(1-pe) 



2 1 -X 



19 



+gauge terms 

Cf i-gn 

+ ln(l - xqx))] 
+gauge terms + 0{e), 



J + +px- (x + 2^—-^^ (Inx - 21n(l - x) 



(99) 



^ _y9a9i _ ^ 

' 27rri 



d(f)dipole „ 

2piPa 

V''{x,xo,e) 

1 

1 - x)-^%l - xoxf J dww-^-\l - w) 





1-e) \ml-Py 



Ca 



-r ( - 1 + X - x^) + 2(1 - pe)^-^'\^s^^ 

VI — X + M / X ZpiQ 



Li -X 

+gauge terms 



,Fi(l,-£;l-2£;- ^) 

V 1 — Xox/ 



— 1 + a; — + 



1 — a; 



X 



+ 



1 1 



e 1 — X 



- (lna; + ln(2-xox)) ^^^3 



X 



1-x + x'^ - 



J- 

1 — a;\ 



ln(l 



— a; I 



1 — X 



+5(1 



^ + In X - 2 ln(l -x)+ ln(l - Xqx)^ ^ 

'1 TT^ 1 

2? ^ 12 ^ 2^ ln(2 - xo) + ln(l - xq) ln(2 - xq) + Li2(xo - 1 



+ 



— x) 



— - ln^(2 — xoj'j + gauge terms + 0{e). 



(100) 



Here 715(5) = 2 denotes the polarizations of a fermion and ns{g) — 2(1 — pe) denotes the 
polariztions of a gluon. Note that the dependece on the momentum q in Vj^"^* and V|"^' 
dropped out after integration. 

6.2 Massive final emitter, massless initial spectator 

Wc now consider the case of a massive emitter (particle k) in the final state and a massless 
spectator (particle a) in the initial state. It is sufficient to consider the case where a heavy 
quark (or antiquark) emits a gluon. The variables z and x of Catani and Seymour are 
given by 



z — 



2pkPa 



2piPa + 2pkPa 

The variables z and w are related by 



X — 



2piPa + 2pkPa - 2piPk 
2piPa + 2pkPa 



(101) 



1 - W, Xo = 1 - 

1 — XqX 



mt 



ml - {pk+Pi-PaY 



(102) 
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Since there is no coUinear singularity, we just have to match the part of the soft singularity 
which corresponds to this dipole factor: 



2(1 -xo) 



X 



1 — z + [1 — X) 1 — X 

We use the same mapping for the momenta as in the previous section: 



(103) 



Pa = XPa, 

Pk = Pk + Pi - (l - X)pa. 



This mapping fulfills the asymptotic behaviour 



lim x = l, \im pk=Pk+ Pi, lira pa = Pa- 



Integration yields 



^dipole 



a.. 



47r/i 



2 \ e 



2piPk ' "'^^ 27r r(l - e) \ml - ^ 



VQg{x,xo,e), 



VQg{x,xo,e) 



a;^(l-x) ^(I-Xqx)" J dww "{l-wY 





C 



l-z + {l~x) 

CpSss'X'^i^ — X)^^'^{1 — XqX 

( 2 



2(1 -xo 



X 



1 — X 

r(2 - 2e) 
-1 



2FAlA-e-2-2e-- 

1 — X V 1 — XqX 



-2(l-Xo) 



X 



1 — X 



(104) 



(105) 



Cf5ss' 5(1 - x) f - (1 + ln(l - xo) - ln(2 - xq)) + 2 - y + ln(l - xq) 
+ ]- ln2(l - Xo) + ^ ln2(2 - Xq) - 2 ln(l - xq) ln(2 - Xq) - 2 Li2(xo - 1) 



+2 ^ln(2 — Xqx) — ln(l — Xqx) — 

+0(5). 



1 



(1 - xo)x^ \ 

1 — xqx y I 1 — X 



(106) 



7 Conclusions 

In this paper we have extended the dipole formalism to processes involving heavy fermions. 
We gave the explicit subtraction terms, together with a mapping of the momenta from the 
{n + l)-parton configuration to the n-parton configuration. We evaluated the integrals of 
the subtraction terms over the dipole phase space to order 0{e^). These ingredients are 
sufficient to set up numerical NLO programs based on the dipole formalism for processes 



involving heavy fermions. An application to single-top production is in preparation |30 



Furthermore we used a systematic approach for the integration and subsequent expansion 
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in e. This method allows us in principle to obtain higher orders in e. Our results are 
naturally expressed in terms of harmonic polylogarithms. We are therefore confident that 
the techniques developed in this paper will prove useful also in the cancellation of infrared 
singularities at next-to- next-to-leading order (NNLO). 
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